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17 Abstract. The evolution of the surfaces of bodies unprotected by either strong magnetic fields or 

18 thick atmospheres in the Solar System is caused by various processes, induced by photons, 

19 energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the 

20 solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, 

21 with implications for their evolution. Ion precipitation produces neutral atom releases into the 

22 exosphere through ion sputtering, with velocity distribution extending well above the particle 

23 escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms 

24 (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) 

25 with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide 

26 mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This 

27 paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA 



28 and suggests the need for quantitative results from laboratory simulations and molecular physic 

29 modeling in order to understand SHEA data from planetary missions. In the Appendix, referenced 

30 computer simulations using existing sputtering data are reviewed. 

31 1. Introduction 

32 Studying the evolution of the surfaces and atmospheres of bodies in the Solar System is 

33 fundamental to our understanding of the present composition of planetary surfaces and 

34 atmospheres. This endeavor entails finding how the rates of the on-going processes vary as a 

35 function of the space environment. Aside from occasional catastrophic events, such as volcanic 

36 eruptions in a few bodies or occasional collisions with comets and asteroids, surface and 

37 atmospheric changes are caused predominantly by the continuous bombardment of the bodies by 

38 photons, energetic ions and micrometeoroids. Yet the actual effects of these incident fluxes on the 

39 present state of planetary bodies are not well described. To investigate this complex topic, we 

40 propose to begin with a much simpler quest by focusing on the subset of planets, moons and small 

41 bodies that are not protected by either strong magnetic fields or thick atmospheres. For surfaces of 

42 exposed bodies (EB) such as Mercury, Moon, and asteroids, directly exposed to the solar wind, the 

43 alteration of the solid surface and the production of the surface-bound exospheres by the impacts of 

44 the time-varying solar wind (SW) over the last 4.54 Gy constitute a relevant component of space 

45 weathering. For other EB, such as Callisto, Europa and Ganymede of Jupiter, energetic 

46 magnetospheric (MS) ions play the major role in altering the respective surfaces and atmospheres. 

47 Hence, we shall focus on the process of ion sputtering (IS) on EB, i.e. on bombardment by either 

48 SW or MS ions. 

49 In the past, the nature of space weathering has been reviewed in detail and the sputter alteration of 

50 regoliths of outer solar system bodies has been discussed (Hapke, 1986; 2001). Moreover, a 

51 mathematical theory describing the optical effects of space weathering has been derived and applied 

52 to the regoliths of the Moon, Mercury and an S asteroid (Hapke, 2001). Whereas Hapke (2001) 

53 discussed the spectral effects and the melting of minerals caused by space weathering, in this study 



54 we consider as its main specific signature the flux of energetic atoms ejected upon impact of 

55 energetic particles on the surfaces. 

56 Although there are other surface-ejected atoms and molecules, such as those released by thermal 

57 desorption (TD), photon-stimulated desorption (PSD) and micrometeoroid impact vaporization 

58 (MIV), we shall show that IS ejecta produced by the incident SW or MS ions provide a unique 

59 window to observe space weathering of EB. These ions may be partly neutralized and back- 

60 scattered from the surface to space (up to 20% for light ions like the SW major components, see 

61 McComas et al., 2009, and Wieser et al., 2009), but a significant fraction of the incident ions, 

62 increasing with ions atomic mass number, can be implanted on the EB surface while ejecting a 

63 surface atom or molecule. Sputtering products from impacts of keV ions can have energies, peaking 

64 at few eV with a high-energy non-Maxwellian tail, up to at least several tens eV for a refractory 

65 material (Goelhich et al. 2000). We refer to this component of the surface ejecta as sputtered high- 

66 energy atoms (SHEA). At these energies, SHEA emitted from regolith can easily escape the local 

67 gravity (e.g., 0.09 eV/nucleon for Mercury and 0.03 eV/nucleon for the Moon) and be distinguished 

68 from other surface-released products from TD, PSD and MIV, all typically <1 eV. Plainaki et al. 

69 (2010) show that even in the case of icy moons the flux of escaping IS ejectia is significantly higher 

70 than the other products (see also Figure A5 in the Appendix). The energy spectra of SHEA, of 

71 course, strongly depend on the incident flux and surface composition. Being electrically neutral and 

72 energetic, SHEA can escape both the magnetic and gravitational field present between their places 

73 of birth (where sputtering occurs) and a SHEA analyzer onboard either an orbiter or a fly-by probe. 

74 If, on the same spacecraft, the SW or MS ions are monitored by a plasma analyzer and the surface 

75 composition of the exposed bodies (EB) analyzed by IR, X-ray, y- ray or neutron spectrometers, then 

76 the detection of mass and energy distributions of SHEA would provide the missing piece in 

77 determining the magnitude and rate of space weathering of the given EB’s surface as well as the 

78 composition of its surface-bound exosphere. 



79 Recent observations of heavy pickup ions at Mercury by the MESSENGER spacecraft (Zurbuchen 

80 et a/., 2008) and the pickup ions from reflected SW protons at the Moon by the Kaguya spacecraft 

81 (Saito et al. , 2008) have shown most clearly that SW-ion and EB-surface interactions are a link 

82 between the physics of space plasma and of surface-bound exospheres. Furthermore, Chandrayaan- 

83 1 Energetic Neutrals Analyzer (CENA) was, in principle, able to measure neutral atoms of 10 eV to 

84 3 keV (Bhardwaj et al., 2005). This sensor observed an energetic neutral signal from the Moon 

85 surface, interpreted as the product of neutralization and back-scattering of the solar wind, probably 

86 prevailing on sputtering signal at the Moon (Wieser et al. 2009). The results and sensitivity of 

87 CENA could provide an indication for estimating an upper limit of the flux of SHEA around the 

88 Moon. These recent results come, however, from in-situ measurements of the already processed 

89 surface releases. Were remote sensing of the surface via SHEA from the vicinity of Mercuiy and of 

90 Moon with appropriate instrumentation available, then more direct and detailed investigations could 

91 be done on the nature of the surface-plasma interaction under different physical conditions, e.g ., SW 

92 condition, solar radiation effect, magnetospheric condition, and surface property. Such 

93 investigations may be carried-on by monitoring SHEA flux intensity, emitting area extension and 

94 particle relative abundances. The comparison between the ground-based observations and 

95 spacecraft measurements and between pickup-ion and SHEA measurements would resolve many 

96 outstanding issues such as the interplay between ion-sputtering and photo-desorption by solar UV 

97 photons, the relative importance of thermal desorption and meteoroid impact as source mechanisms 

98 of the sputtered exospheric atoms. 

99 Clearly, to accurately interpret any SHEA data from space weathering effects on surfaces of EB will 

100 require active participation of physicists doing sputtering experiments in laboratories directly 

101 applicable to the interactions between SW or MS ions and EB surfaces. Only such experiments can 

102 quantify the microscopic processes controlling the sputtering yield Y*, the number of released 

103 particles per incident ion, basic to remote sensing in SHEA. 



104 Remote sensing EB via SHEA by orbiters or fly-by probes can also provide information to 

105 complement the observations from Earth or by instruments landed on these solid bodies. Although 

106 fly-by missions offer only brief observation of one body, each mission could be planned to fly by 

107 several bodies. The advantage of orbiters over landers, besides cost, is its global survey under 

108 varying conditions over longer time periods. In the case of orbiting larger planets with many 

109 moons, the ability to observe several moons has been successfully demonstrated by missions 

110 Galileo and Cassini. These and other orbiter missions, unfortunately, are not equipped to study 

1 1 1 space weathering of the EB. To examine the issue of SHEA capability on future EB missions, an 

112 in-depth discussion is necessary. 

113 To begin this discussion, we start with the data and techniques currently available to assess whether 

1 14 or not SHEA instruments are critical to future orbiter or fly-by missions. To this end, details are 

115 presented in the following manner: the production of SHEA in §2, justification for SHEA 

116 observation in §3, the need for laboratory-based ground truth in sputtering in §4, and the 

117 conclusions in §5. Examples of computer simulations of SHEA emissions from Mercury, Moon, 

118 asteroids and Jovian Moons, based on existing data and theories, are presented in the Appendix. 

119 

120 2. Production of SHEA 

121 The IS results from the impact of an ion of mass m\ onto a solid surface. If the ion incident energy 

122 Ei is high enough, surface atoms may be ejected. Some IS processes producing SHEA are 

123 represented in Figure 1. For oblique incidence, ion sputtering can be a single-step process, often 

124 called “knock-on”, in which the ions directly eject surface atoms (a). Otherwise, a multi-step 

125 process takes place, often called “collision cascade” (d). Light incident ions are often backscattered 

126 in layers near the surface, and occasionally they may be neutralized in the process before returning 

127 to space - not shown, but would be like in (e) without the second collision. Backscattered ions can 

128 trigger a cascade of collisions among atoms close to the surface. While the heavy incident ions 

129 produce forward-directed recoils. 



130 The energy transferred in the first collision to a surface atom is given by classical mechanics: 


T = T m coia r 
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132 where E, and m\ are the incident ion energy and mass, respectively, m 2 is the mass of the struck 

133 atom (the recoil), T the energy transfer, T m its maximum value, and <Xr the scattering angle of the 

134 recoil atom (see Figure 1). Collisions below the surface layer involve both the projectile and the 

135 recoil atoms, with the cascade of collisions eventually leading to sputtering, i.e. the ejection of an 

136 atom or molecule from the solid. For a regolith material (independently of composition or porosity), 

137 the ejected particles are mostly neutral atoms (Hofer, 1991). For ejected atoms or molecules of 

138 species n with partial sputtering yield Y„, the normalized distribution of ejecta (/$„) from a 

139 refractory material, as a function of ejecta energy E e , can peak at few eV (Gnaser 2007; Hofer 1991) 

140 and can often be empirically reproduced by the following function (Sigmund, 1969; Sieveka and 

141 Johnson, 1984): 

cos«„ (2) 

143 where Eb, n is the surface binding energy of the ejected atoms, a n the polar angle of the SHEA with 

144 respect to the surface normal (Figure 1), and c„ the normalization constant. Gnaser (2007) showed 

145 that the effective binding energy Eb,„ is typically lower than the bulk cohesive energy. For 

146 refractory materials, the difference between the two can be as much as 50%, but more typically ~ 

147 10-20%. For volatile materials that dominate the outer solar system, the difference can be an order 

148 of magnitude (e.g., Reimann et al. 1984; Johnson et al. 2010). Empirically, all variables in Eq. (2) 

149 can be measured, except Eb,„. By fitting Eq. (2) to laboratory data on sputtering, therefore, can 

150 uniquely determine Eb, n . 

151 Samples for comparing computed with measured _/s> as functions of E e for different incident ions 

152 and solid targets are shown in Figure 2. Panels (a) and (b) are fs,„ computed for Na ejected by 




153 protons incident on a planetary-like mineral for different values of Ebjfa and of Et, respectively, 

154 using Eq. (2) averaged over angle a„ (Sigmund, 1969; Sieveka and Johnson, 1984). It is clear that 

155 Ej sets the upper limit on E e , while £* affects the energy at which the distribution peaks. Panel (c) 

156 shows good agreement between Monte Carlo SRIM (Ziegler et al., 1966) simulation results (for a 

157 surface composition derived by Goettel, 1988) and Eq. (2) in the high-energy tail in the case of 1- 

158 keV protons on a planetary-like surface. Moreover, panel (d) compares Eq. (2) with experimental 

159 results of Ar + impacting on W at four different values of E> for 9, - «„ = 0 (Figure 3 of Goelhich et 

160 al., 2000) with Eq. (2); the agreement improves for is, > 500 eV. The spectrum of the ejected Na 

161 shown in panel (e) is converted from velocity to energy E e as the independent variable, resulting 

162 from bombarding a Na 2 SC >4 target with 3.5-keV Ar~, as might be the case for surfaces of Io (Wiens 

163 et al. 1997) or, possibly, certain regions of Europa although Na is often in an ice matrix (e.g., 

1 64 Johnson et al. 2002). The spectrum fits the form of Eq. (2), which has a measured tail extending to a 

165 few eV, but peaks at ~ 0.3 eV, well below that shown in panel (d). They also showed that the 

166 Monte Carlo SRIM is able as well to reproduce the process for different impact energies and angles. 

167 Panel (f) gives the energy spectra of sputtered D 2 O and SO 2 from 5-keV Ar + impacting a heavy- 

168 water ice matrix containing SO 2 (Johnson et al., 2010). Figure 2 demonstrates the wide 

169 applicability of Eq. (2), except for the lowest energy portion showed in panel (f), as explained in 

170 Johnson et al. (2010), and the need to establish £*,„ for incident ions and targets relevant to die 

171 study of space weathering of selected EB. 

172 The angular distribution of the ejected atoms depends on the incident ion mass so that a general 

173 expression is not easily defined; detailed discussions can be found in Hofer (1991) and Gnaser 

174 (2007). For heavy incident ions, the ion impact direction does not have a large effect on the 

175 distribution in ejection angle a„, which is often approximated cos*(a„), where k is usually between 

176 1, as in Eq. (2), and 2. For light ions, the angular distribution is related to the ion impact direction, 

177 and exhibits a maximum close to the mirroring angle. For a surface composed of a number of 

1 78 different atomic species, the angle-averaged differential flux of sputtered atoms is: 
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where C„ is the relative surface abundance of the atomic species n, and 4>/ is the incident ion flux. 
The total sputtering yield Y=\(d&/'dEe)dE e , in general, depends on the impinging ion mass and 
energy and on the surface mineralogy. Averaged over the solar-wind ion energies, Y can range from 
0.01 - 0.1 (Lammer et al., 2003; Johnson and Baragiola 1991) for refractory surfaces, whereas it 
ranges in between 10 - 1000 for icy surfaces of the Jovian moons when bombarded by hundreds 
keV MS heavy ions (Johnson 1990; Fama et al. 2008; Johnson et al. 2010). These values are 
reduced by the regolith porosity (Cassidy and Johnson, 2005). The yield is also a function of the 
incident ion’s mass and nuclear charge. In general, every precipitating ion contributes to sputtering 
from the EB surface. For example, accounting for the solar wind abundance of the ions, the net 
sputtering rate generated by protons with respect to other solar-wind components, like a particles or 
high-charge-state particles, is estimated to be comparable (Johnson and Baragiola, 1991). In the 
case of the icy moons of the giant planets hit by heavy and energetic magnetospheric plasma ions, 
the ejecta are dominated by low-energy atoms and molecules. Since the yields from such surfaces 
can be large, both simulations (Cassidy et al. 2009) and experiments (Johnson et al. 2010) show that 
the trace species are carried off with the ice matrix. The two sets of EB, one exposed to SW only 
and the other exposed predominantly to magnetospheric plasma, make a comparative study that 
would improve our understanding of the mechanism of space weathering. 


3. The uniqueness of SHEA observation 

3. 1 Selecting a starting point in the study of surface evolution 

The understanding of the role played by SW and MS ions, solar radiation and micro-meteorites in 
bombarding, in altering the surfaces and atmospheres of these bodies, as well as the determination 
of the mass loss rate of the respective bodies (Killen and Ip, 1999; Madey et al. 2002) provides a 
relevant contribution to the study of the evolution of solid bodies of the Solar System. To begin this 



204 ambitious and challenging study, we have, as stated in §1, selected the EB in the Solar System, that 

205 are not protected by either strong magnetic fields or thick atmospheres. Such bodies are directly 

206 exposed to the incident radiations, and the resulting released atoms and molecules can escape with 

207 least hindrance. On the other hand, those atoms that fail to escape populate the surface-bound 

208 exospheres ( e.g Johnson, 2002). The choice of EB also minimizes interference, such as deflection 

209 by strong local magnetic fields or/and scattering by intervening atmospheric particles, on the 

210 incident radiation as well as on the ejecta from the site of impact. 

21 1 Among the processes occurring on the surfaces of EB, which include TD, PSD, IS and MIV, we 

212 select, also stated in §1, IS (ion sputtering) the first process for investigation. Our choice of IS on 

213 EB to begin our study on surface evolution is not just because we recognize the principal role of the 

214 time-varying ion flux intensity over the last 4.54 Gy in space weathering of bodies in our Solar 

215 System (Orsini et al, 2009c), but also due to the fact that three necessary sets of observables can be 

216 made accessible. These complementary observables are: the incident radiation, which has been and 

217 will continue to be monitored by planetary missions; the surface composition and mineralogy of 

218 EB, which have been and should be investigated by space-borne X-Ray, IR, neutrons, and gamma- 

219 ray spectrometers; and the ejecta of IS, which have distinct features that favor direct and precise 

220 detection and analysis, but yet to be implemented. Recently, Kaguya and Chandrayaan-1 spacecraft 

221 had X-Ray, IR, neutrons, gamma-ray and particle analyzers. Although Kaguya particle analyzers 

222 measure only ions and electrons, some of SHEA are ionized. A joint analysis of these kinds of 

223 observations could provide hints in this study. The future BepiColombo Mission, already including 

224 in its payload all these sensors and especially a dedicated SHEA detector, promises outstanding 

225 outcomes (see Section A 1 ). 

226 3.2 SHEA detection for observing space weathering 

227 a) Energetically distinct 

228 Different release processes produce particles within different energy ranges (Wurz and Lammer, 

229 2003; Milillo et al., 2005; Leblanc and Johnson 2010). The ejected atoms and molecules, depending 



230 on their velocity, can either return to the surface, become part of the atmosphere, escape the 

231 gravitational field, or be photo-ionized and picked-up by planetary magnetospheres. The velocity 

232 distributions are different for the relevant processes, and thus can serve as important signatures of 

233 the processes involved. TD and PSD are more effective for volatiles (like H, He, Na, K, S, Ar) and 

234 have typical energy below 1 eV (dashed lines in Figure 3(a), (b) refer to 2-eV-Na, that is, the escape 

235 energy at Mercury), while IS and MIV are effective also for refractory species (e.g.. Mg, Al, Si, 

236 Ca), thus producing more energetic ejecta closer to stoichiometric composition. In contrast to the 

237 MlV-released particles having a Maxwellian distribution of an expected peak corresponds to 

238 -2500-5000 K (Eichhom, 1978) or a peak particle energy of -0.6 eV, the high-energy tail of IS 

239 ejecta, SHEA, on the other hand, can in principal have surface release energies above 10 eV 

240 (Gnaser 2007; Wiens et al. 1997), more than sufficient to escape the local gravity (e.g., 0.09 

241 eV/nucleon for Mercury, 0.03 eV/nucleon for Moon). This means that releases from all other 

242 processes can be excluded, when analyzing IS products through SHEA detection (Figure 3(c) and 

243 also Figure 2). Nevertheless, the escape fraction of released particles depends on each specific case 

244 (escape velocities, main release processes, surface properties, external conditions) and it is a 

245 complicated quantity to estimate. Generally, one of the main processes responsible for the total 

246 surface material loss rate is IS, but minor contribution can be due to radiation pressure for specific 

247 species and to the other release processes, as well. 

248 SHEA may also be distinguished from back-scattered atoms (BSA). This population is not a 

249 negligible fraction of material leaving the surface, but definitely not of surface composition. BSA 

250 are just neutralized impacting ions that are reflected back from the surface, so that their energy is 

251 comparable to that of the incident ions. Back scattering is much more efficient for light species, like 

252 H, so that both their flow velocity and energy are well separated from those characterizing the IS 

253 ejecta, SHEA. This means that an instrument is able to discriminate between these two signals 

254 provided that its ToF or energy resolution is high enough. If we consider 1 keV proton onto regolith 



255 surface, we can assume a total yield of ion sputtering about 10%, and that of backscattering about 

256 20%, then, the expected fluxes are comparable. 

257 SHEA mapping on EB is distinctly different from ENA (energetic neutral atoms) imaging remote 

258 plasma such as planetary magnetospheres or moons (e.g., Hsieh & Curtis, 1988; Krimigis et al., 

259 2004). The latter relies on the production of energetic atoms by charge-exchange between energetic 

260 ions and ambient atoms and molecules along the line of sight and within the solid angle of the ENA 

261 imager. The intensity of charge-exchange ENA flux is, therefore, a column-density measurement 

262 along the line of sight. The choice of EB as the solid target and IS as the process effectively renders 

263 any ENA produced along the line-of-sight between EB and the SHEA detector insignificant, 

264 because the charge-exchange cross section typically < 10~ 14 cm 2 for ~10-eV ions (e.g., Lindsay and 

265 Stebbings, 2005), the number density of atoms in interplanetary space is ~10 1 cm -3 (e.g., Bzowski, 

266 Fahr & Rucinski, 1996, Bzowski et al., 2008), and distance between EB surface and the observing 

267 spacecraft (s/c), hence the path length for ENA, ~10 2 ' 3 km. Hence, the product of these three 

268 quantities indicates that the maximum ENA flux that can reach the observing s/c from the observed 

269 EB would be -10 -6 of the ambient ion flux. This is orders of magnitude smaller than the expected 

270 SHEA flux under the bombardment of the same ion flux, due to the fact that the all-species- 

271 integrated sputtering yield is of the order 0.1, in the case of a regolith surface hit by 1-keV proton. 

272 Moreover, ion fluxes at 10 eV are usually negligible in the EB environment; generally, charge- 

273 exchange ENA are in the keV range, when the plasma is mainly SW, or they can have higher 

274 energies when considering the giant planets magnetospheres, and ion directions are generally not 

275 from the body to the s/c. So the expected ENA flux comes from different directions and at different 

276 energy range from those of SHEA. 

277 Having distinguished SHEA from back-scattered neutrals, charge-exchange ENA and ejecta of 

278 other surface-altering processes, we arrive at the unique advantage of observing targeted EB in 

279 SHEA. 


280 


b) SHEA-mapping: instantaneous & localized 



281 While the ground-hugging exospheres of EB maybe a mixture of lingering releases from all other 

282 surface processes over time, escaping SHEA, on the other hand, travel ballistic trajectories from 

283 their ejection sites or ion-impact site to the observing spacecraft, thus carrying instantaneous and 

284 localized information on their origins. SHEA enable us to directly map the spatial distribution of the 

285 ion impact flux in time. Correlating observed time profile of SHEA with that of the impinging ions, 

286 e.g . SW or MS ions, bombarding the surface, with the knowledge of surface composition provided 

287 by means mentioned in the beginning of §2, it is not difficult to imagine how the specific yield and 

288 erosion rate could be obtained, within the time-spatial and mass resolution of the SHEA instrument. 

289 4. Necessary ground truth 

290 In the face of the attractive and unique advantage of observing surface erosion of EB by IS via 

291 SHEA, we caution the need for minimizing the uncertainties from the complexity of the surface 

292 being bombarded by ions of different species and energy and ejecting SHEA of different species 

293 and energies. This prerequisite for extracting information reliably from the three sets of data - 

294 incident ion fluxes, surface composition, and SHEA maps - must be guided by solid ground truths 

295 found only in extensive laboratory data on sputtering mechanisms and yields. 

296 Quantitative laboratory simulations and computer modeling of IS occurring on EB are essential for 

297 understanding SHEA data from planetary missions. This is analogous to the need for ground truth 

298 in remote sensing: only on-site measurements that help calibrate aerial photographs and satellite 

299 imagery can make data interpretation and analysis credible. It has been suggested that 

300 measurements of composition and kinetics of atoms and molecules in a body’s exosphere during an 

301 orbiting mission could determine the importance of the different surface release mechanisms, and 

302 the surface composition (e.g., Johnson et al., 1998). In particular, the IS process will require 

303 laboratory measurements to support existent and future planetary missions. We cite here some 

304 existing use of theoretical knowledge of IS on space data (see Appendix), and where laboratory data 

305 are needed to resolve complexities that theory alone proves difficult. 



306 As mentioned in §2, the ejection of surface atoms or molecules by IS is characterized by the yield Y. 

307 At projectile energies of the order keV, e.g., SW, IS occurs due to electronically-elastic knock-on 

308 (ballistic) processes that are fairly well described by the linear cascade theory (Sigmund 1969). For 

309 certain insulators, the electronic excitations produced by the projectile can live long enough to 

310 produce what is known as electronic sputtering (Johnson, 1990). The relative importance of these 

3 1 1 two processes depends on ion velocity and charge state. 

312 According to the standard linear collision cascade theory, the elastic sputtering yield for atomic 

313 targets is proportional to the ratio of the energy deposited at the surface and the binding energy of 

314 the surface atoms. This concept has been widely used to estimate the contribution of sputtering to 

315 the exosphere of Mercury, the Moon, the NEO and the Jupiter moons (Wurz and Lammer 2003; 

316 Massetti et al., 2003; Leblanc et al. 2007; Mura et al., 2009, Plainaki et al., 2009, Plainaki et al., 

317 2010). Concerning this last case, the sputtering is much more complex, since electronic sputtering 

318 occurs. In fact, this process is very effective in materials with low cohesive energies such as the 

319 frozen gases in the outer solar system, often referred to as “ices”. For such materials the total 

320 electronic sputtering yield Y is often proportional to the square of the electronic-stopping cross 

321 section. Early laboratory sputtering data by Brown et al. (1982) were used to predict the principal 

322 atmospheric component of Europa, O 2 , and its average column density (Johnson et al., 1982). In 

323 addition, the predicted large sputtering yields have led to the suggestion that other trace species 

324 should be present (Johnson et al. 1998). 

325 Electronic sputtering is closely related to desorption induced by electronic transitions (DIET) (e.g., 

326 Madey et al. 2002). In DIET, an incident electron, ion or photon excites a surface state, which can 

327 relax by ejecting an ion or a neutral. This is a process that is linear in the excitation cross sections 

328 and is responsible for the sodium atmospheres on Mercury and the Moon (Y akshinskii and Madey 

329 2000). The DIET process occurring on Mercury and the Moon is molecularly specific and is, 

330 therefore, efficient for specific trace species (primarily the alakali’s) or molecules adsorbed on 

331 refractory surfaces. Energetic electrons, ions or photons can produce deep excitations which, in 



332 principal, can result in the ejection of a large variety of surface species. However, such excitation 

333 events typically occur with lower probability. The interest in knock-on sputtering is that it is more 

334 robust and could eject into the gas phase species that are more representative of the surface 

335 composition. With the discovery of calcium ejected from Mercury’s surface, this would appear to 

336 be bom out. However, predictions for bodies with silicate surfaces, like the Moon (Johnson and 

337 Baragiola 1991), have been much less successful than is the case for the icy bodies discussed above. 

338 This fact is primarily because the yields are small («1 atom per ion), e.g., for SW bombardment, 

339 so that the sputtering of an element is more sensitive to its molecular surroundings, and, as 

340 discussed below, there is insufficient data on refractory planetary materials. In attempting to model 

341 this process, there are several reasons that would discourage the use of the linear cascade theory to 

342 estimate the elastic sputtering yield contribution to the planetary exospheres as currently being 

343 applied in atmospheric models. The theory was developed for mono-atomic targets, it assumes a 

344 constant binding energy for atoms at the surface, and since it is based on a transport theory 

345 approximation, it only works for amorphous materials. Of critical importance in planetary science is 

346 the so-called “threshold regime”, where the model breaks down and empirical models are used. 

347 Sputtering becomes even more complex, if the target consists of two or more different atomic 

348 species. The complication arises because the energy transfer from the projectile to the various target 

349 species is different. More important, each species has a different binding energy to the lattice and, 

350 therefore, irradiation leads to enhanced diffusion and depletion of the more volatile species resulting 

351 in a change of the composition of the solid with depth. 

352 Sputtering yields are usually measured on relatively flat laboratory surfaces. However, meteoritic 

353 bombardment over millions or billions of years on the surface of an astronomical body produces a 

354 regolith, a porous surface composed of grains formed by cumulative fracture and crater ejecta. Ions 

355 impact a regolith structure over a range of incident angles. Since the sputtering yield depends 

356 steeply on the local incidence angle 0 (for ices the standard linear cascade theory predicts a 

357 dependence of cos' f 0, where f is between 1-2 and is nearly independent of the projectile energy; 



358 Fama et al., 2008), one would expect that the yield from a regolith would be different compared to a 

359 hypothetical flat surface. This effect has been evaluated using Monte Carlo simulations by Cassidy 

360 and Johnson (2005), who found that the total sputtering could be significantly lower than the 

361 laboratory yields. 

362 In contrast with the numerous studies of the sputtering of water ice (see Baragiola et al., 2003 and 

363 Fam£ et al., 2008 and references therein) which have been of useful application for analysis of outer 

364 planetary systems (e.g., Johnson et al. 2008), there are few laboratory measurements of sputtering 

365 of neutrals and secondary ions from minerals (Betz and Wehner, 1983; Jull et al. 1980; Elphic et al. 

366 1991; Betz and Wien, 1994; Wiens et al. 1997). Therefore, measuring sputtering rates and velocity 

367 distributions of sputtered species from minerals and ices relevant to planetary surfaces is essential to 

368 support SHEA data from future planetary missions. Because such measurements are time intensive 

369 and can often not be made over the full energy range required, simulations of sputtering are critical 

370 for extending the range of applicability, especially in the threshold region. Both Monte Carlo test 

371 particle simulations and Molecular Dynamics simulations have been carried out. The Monte Carlo 

372 simulations, typically only track recoils with energies much greater than the cohesive energy of the 

373 solid, and necessarily give results equivalent to those obtained from the linear Boltzmann equations. 

374 The best known of such calculations are the heavily used TRIM/SRIM models (see Section 2). 

375 However, these are applicable only in regions in which linear cascade model is valid and fail in the 

376 threshold regime. Much more useful are the Molecular Dynamics methods which are, of course, 

377 much more computationally intensive. In such models the atoms and/or molecules in the materials 

378 interact with each other and with the incoming ions via intermolecular potentials. To date, they have 

379 been primarily applied to model materials (Tucker et al. 2005; Bringa et al. 2000), but extensions to 

380 materials with compositions and properties similar to surfaces of the Moon and Mercury are 

381 feasible. Because the intermolecular potentials are not known in details for complex materials, both 

382 types of simulations are typically calibrated to experiment. Therefore, they are primarily useful as 

383 means for extending the range of the available data. This combination of laboratory experiments 



384 and numerical simulations will be essential to support the proposed SHEA instrumentation and 

385 mission design. 

386 

387 5. Conclusion 

388 In the interest of understanding what kind of evolution led to the present composition of planetary 

389 atmospheres, including that of Earth, we need to learn how the current rates of the on-going 

390 processes that cause surface modifications and particle escapes vary under different conditions. To 

391 begin, we suggest concentrating on planets, moons and small bodies that are not protected by either 

392 strong magnetic fields or thick atmospheres, i.e. EB (exposed bodies). Furthermore, we identified 

393 IS (ion sputtering) on EB as the first process for detailed investigation for three reasons: 

394 1. Incessant bombardment by either SW (solar-wind) or MS (magnetospheric) ions on the 

395 respective EB constitutes predominant relevant process in altering the surface and consequently the 

396 associated ground-hugging exosphere. 

397 2. The ejected products of IS on the surface, SHEA (sputtered high-energy atoms), are mostly 

398 neutral and energetic enough not only to be distinct from surface-released particles from other 

399 processes, but also to escape local gravity and magnetic field for remote sensing. 

400 3. Remote sensing in SHEA can provide mapping of the EB trader IS with temporal and mass 

401 resolution. 

402 We illustrated what is possible and what is needed to realize remote sensing IS on EB in SHEA. In 

403 the Appendix, referenced computer simulations using existing sputtering data and realistic EB 

404 models are shown. We hope this paper has initiated the drive towards determining how the rates of 

405 the on-going IS processes that cause the changes vary under different SW and MS conditions. 

406 For such an effort, we see that parallel to planning flyby or orbiter missions to EB and developing 

407 SHEA instrumentation (e.g.: BepiColombo/SERENA/ELENA Orsini et al. 2009, 2010) for such 

408 missions, performing sputtering experiments in the laboratories using appropriate incident ions and 

409 impacting surface analogues is equally necessary. The latter would indeed produce data crucial for 



410 the planning of the missions and design of SHEA instruments as well as for extracting factual 

411 information from the ensuing SHEA data. 

412 
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417 Appendix: Simulations 

418 Based on currently available IS data and theoretical models of EB, we present here material 

419 extracted from recently published papers, on what could be expected from SHEA imaging, by 

420 simulation of the following EB in their particular environments: A. I) Mercury, when the 

421 interplanetary magnetic field (IMF) configurations and SW conditions permit the plasma to reach 

422 the planetary surface, A. 2) the Moon, in the SW while outside the Earth magnetosphere, or when it 

423 crosses the plasma sheet, A. 3) asteroids and small bodies continuously exposed to the solar-wind 

424 plasma, and A. 4) the Jovian moons, Europa and Ganymede, when embedded in the Jupiter 

425 magnetosphere radiation belts. 

426 A.l Mercury 

All SHEA measurements of Mercury should be particularly intriguing, since they would give the 

428 opportunity to investigate the MS and planet interaction with the intense SW flux at about 0.3 AU. 

429 Many authors (e.g., Killen et al. 2001; Sarantos et al. 2001; Kallio and Janhunen, 2003; Massetti et 

430 al., 2003) showed that under specific IMF configurations, SW can enter through the cusps in the 

431 dayside MS, eventually reaching the surface at mid-latitudes. Under different IMF conditions, the 

432 configuration of the Hermean MS changes so that the area of high proton precipitation (hence: of 

433 subsequent SHEA release) moves accordingly (e.g., Sarantos et al., 2001; Kallio and Janhunen, 

434 2004; Massetti et al., 2007). Ground-based observations, indeed frequently indicate mid-latitude 

435 exospheric Na densities to vary over timescales of hours (Leblanc and Johnson et al. 2010). The 



436 nature of such variations is still unresolved, but it appears to be related to plasma precipitating 

437 regions (Orsini et al, 2008; Mura et al. 2005; Leblanc et al., 2007; Mura et al., 2009). 

438 The main constituents of the Hermean exosphere are probably volatiles released thermally or by 

439 PSD. Not contained in the exosphere, but directly ejected and escaping the gravity will be the flux 

440 of SHEA. The sensor ELENA (Orsini et al. 2009), part of the SERENA particle package (Orsini et 

441 al. 2010) on board BepiColombo ESA-JAXA mission (launch 2014; Benkhoff et al., 2010) will 

442 permit for the first time to map the IS emission, less intense than the PSD emission, but more 

443 effective in releasing refractories from the Hermean surface (Milillo et al. 2005). The flux of -1-5 

444 keV SW protons hitting the Hermean surface is estimated ~10 9 cm" 2 s' 1 ; a total sputtering yield 

445 -10% of the incident ion flux would lead to a total sputtered flux -10 8 cm" 2 s' 1 . Approximately 50% 

446 of the ejected particles escape the planet along ballistic trajectories; and -1% of these particles have 

447 enough energy (> 20 eV) to be detected by the ELENA sensor. For comparison, the back scattering 

448 flux is of the same order of magnitude as that of the sputtered signal, but with an energy spectrum at 

449 higher energies. The top-left panel of Figure Al shows simulated total sputtered flux from 

450 Mercury’s surface over the Northern hemisphere (Mura et al., 2005). The portion of the surface 

451 seen in SHEA from a vantage point at 400 km altitude, latitude 45° and LT 1200, is illustrated in 

452 the lower-left panel. 

453 As BepiColombo Mercury Planetary Orbiter (MPO, where SERENA will be mounted) will fly over 

454 Mercury at low altitudes (orbit: 400x1550 km, see Benkhoff et al., 2010), ELENA’s narrow field of 

455 view (4° x 76°, with 4° x 4° resolution) will ground-frack the Mercury’s surface in SHEA along the 

456 MPO orbital path, as shown in the upper-right panel of Figure Al. The lower right panel puts a 

457 single scan in perspective to the SHEA emitting region shown in the lower left panel. 

458 Eventually, SHEA imaging by ELENA will map the proton precipitation regions, with the help of 

459 the simultaneous detection of the back-scattered neutrals, with surface spatial resolution between 1 5 

460 and 50 km, depending on s/c altitude. Simulations show that ELENA’s spatial and time resolution 

461 capability will allow monitoring the dynamical behaviour of the magnetospheric configuration; 



462 whereas its spatial resolution of tens km will allow to discriminate surface emissivity variations: as 

463 explained in §4, the intensity of the directional SHEA signal depends on both ion precipitation flux 

464 and surface properties, like composition and intrinsic sputtering yield. 

465 A.2. The Moon 

466 The relative importance of IS as a source process for the lunar exosphere remains questionable, 

467 despite the extensive telescopic observations of Na on the Moon. This is due to the limitations of 

468 the viewing geometry from Earth. The Na emission clearly decreases when the Moon lies inside the 

469 terrestrial magnetosphere as would be consistent with a substantial reduction of ion precipitation 

470 contribution to this process (e.g.. Potter and Morgan, 1994), but this decrease coincides with a 

471 changing line of sight: the observer at or close to full Moon is limited to observing the terminator 

472 limb, where flux-dependent sources are weak. As a result, interpretation of the existing ground- 

473 based data must necessarily rely on transport models, which account for this geometrical effect. 

474 Solar-wind proton sputtering of Na (Figure A2) is suggested to be unimportant under the 

475 assumption of yields 0.01 -0.1 per ion, contributing ~ 1 % of the total sodium exosphere. Based on 

476 the yields for desorption induced by electronic transition processes that arc measured in laboratory 

477 studies (e.g., Madey et al., 1998; 2002), PSD is the dominant sodium source process, while 

478 vaporization of regolith material caused by MIV has been suggested to constitute up to 50 % of the 

479 local density of sodium at the terminator, although models may disagree (e.g., Leblanc and Johnson 

480 2010). For refractory species, which remain undetected as neutrals, yet have been identified as 

481 pickup ions (Mall et al., 1998; Yokota et al., 2009), impact vaporization is expected to be the 

482 dominant source process, although ejection in the form of molecular oxides and their subsequent 

483 photo-dissociation remains a candidate (Morgan and Killen, 1997). However, considerable 

484 uncertainties still exist regarding both the flux of micrometeoroids at 1 AU (Cintala, 1992; Love 

485 and Brownlee, 1993) and the effect of multiply-charged solar wind heavy ions, which are common 

486 during Coronal Mass Ejection (CME) events. In-situ SHEA measurements are needed to constrain 

487 the sputtering source. Furthermore, the detection of asymmetries in back-scattered fluxes linked to 



488 magnetic anomalies (observed by Wieser et al., 2010) allows remote sensing of the effects of SW 

489 interaction with the micro-magnetospheres due to local magnetism. Simultaneous observation of IS 

490 emission in SHEA would add necessary details on surface-SW interaction. 

491 The sputtered density of the major species of the lunar regolith (O, Al, Si, Mg, ...) is shown in Fig 

492 A2(b). Sputter yields and elemental abundances consistent with Wurz et al. (2007) were assumed. 

493 At 50-km where a possible orbiter could be located, all the IS species (mainly O and Si) have 

494 densities of the order of 4 cm" 3 in agreement with the estimation of Milillo et al. (2010). This flux 

495 might increase following the exposure of the lunar surface to Earth’s plasma sheet ions, which 

496 impart approximately three times higher energy than SW ions. Measured enhancements of the 

497 sodium exosphere at full Moon have been correlated to such plasma sheet crossings (Wilson et al., 

498 2006; Sarantos et al., 2008). 

499 A3. Asteroids and small bodies 

500 Asteroids suffer erosion and surface alteration from SW and solar and galactic cosmic-ray 

501 bombardment, as well as from solar-photon irradiation and micrometeorites gardening. 

502 Consequently, the relevant surface release processes, when they are within few AU from the Sun, 

503 are IS, PSD and MIV. TD is strongly temperature dependent; hence, its contribution to exosphere 

504 generation becomes important only at about 1 AU from the Sun, and increases when moving 

505 towards perihelion. The detection and analysis of SHEA from asteroids separates IS from the 

506 contribution from other release processes, thus SHEA detection would enable speculating on the 

507 surface erosion under different environmental conditions. SW sputtering investigation provides 

508 important clues on the evolution of a planetary body. 

509 Solar-wind precipitation on the surface of an asteroid can be strongly influenced by the presence of 

510 magnetic fields. This seems to be the case of Vesta; in fact, Vemazza et al. (2006) identified a lower 

511 limit of 3 10 16 A m 2 for eventual possible Vesta magnetic dipole, capable to deviate the solar wind 

512 away from its surface. Not only a dipole can deviate SW from hitting the surface of an asteroid, but 

513 also smaller magnetic structures, known as mini-magnetospheres (Winglee et al., 2000), similarly to 



514 what has been imaged at the Moon in back-scattered neutral atoms, showing a reduction of neutral 

515 flux from the surface corresponding to a strong magnetic anomaly (Wieser et al., 2010). The 

516 possible presence of such magnetic structures can cause a reduction of the SHEA flux released from 

517 an asteroid, thus minimizing local erosion and surface alteration effects. 

518 A study on asteroids’ exosphere based on the simulation of the various release processes on the 

519 surface of the body has been performed by Schlappi et al. (2008), for the asteroids (2867) Steins 

520 and (21) Lutetia, in preparation of the Rosetta flybys. They found SW sputtering to be the most 

521 important exospheric supply process on the sunlit side of an asteroid. At the near Earth distances, IS 

522 is expected to be even more significant. 

523 The escape velocity of a Near Earth Object (NEO) is very low (i.e. 0.52 m/s for a NEO of mass 

524 ~10 12 kg and of radius ~0.5 km), the particles released from the surface of a NEO are, therefore, 

525 essentially lost in space. Given a specific model for the simulation of the various release processes 

526 happening on the surface of a NEO, the efficiency of each of the particle-release processes can be 

527 estimated. Clearly, identifying the NEO surface properties and its interactions with SW can provide 

528 important information on the effects of space weathering on localized surface regions as well as the 

529 global evolution history of the body. 

530 Plainaki et al., (2009) applied the Monte Carlo SPAce Weathering on NEOs (SPAWN) model to 

53 1 obtain the sputtered distribution around a NEO as a result of its exposure to SW (Figure A3). They 

532 found that significant sputtered fluxes could reach a maximum value of 10 n particles m' 2 s' 1 around 

533 the NEO. The major component of sputtered species is expected to be H. The simulated density, 

534 produced by all species of sputtered particles emerging from a NEO surface, is calculated to be ~3 

535 TO 6 particles m' 3 near the NEO surface. The expected SHEA (E > 10 eV) fraction results in ~1% of 

536 the total released particles. On the other hand, the contribution to the total density of the volatiles 

537 emerging from the NEO surface, via the PSD process, is ~ 1-10 8 particles/m 3 



538 The global analysis of the sputtering erosion of the NEO surface would provide unique information 

539 about the present and the past of the NEO’s surface, revealing the mechanism through which the 

540 solar wind has interacted with the surface atoms, in the past millions of years. 

541 A4. Europa and other Jovian moons 

542 The radiation environment of Europa consists of intense H + , 0 + , S', and C + ion fluxes, in the energy 

543 range from keV to MeV (peaking at ~ 100 keV). These ions can erode the surface of Europa via ion 

544 sputtering, ejecting up to 1000 H 2 O molecules per incident ion, and also break the chemical bonds 

545 of the ejected species resulting in the formation of new molecules (e.g. O 2 ), a process called 

546 radiolysis. The neutrals produced have a characteristic spectrum (Cooper et al. 2001; Strazzulla et 

547 al., 2003; Paranicas et al., 2002). Plainaki et al. (2010) found that the most significant sputtered- 

548 H 2 O emerging flux and density come from impinging S + ions, and they amount to 66% and 59% of 

549 the total (3.210 13 H 2 O m' 2 s' 1 and 2.7T0 10 H 2 O /m 3 , respectively). The total sputtering rate for 

550 Europa was calculated to be ~ 10 27 H 2 O/S with escaping ratio 22%. This value, locally on the 

551 moon’s surface, may exhibit variations; probably, it is higher in the trailing face, where the 

552 precipitation is foreseen to be more intense. In fact, this result is inside the range for the Europa loss 

553 rate given in literature and ranging between a few 10 26 H 2 O s' 1 and 10 28 H 2 O s' 1 (Lanzerotti et al., 

554 1982; Johnson et al., 1981; Eviatar et al., 1981, 1985; Shi et al., 1995; Ip, 1996). A similar result is 

555 also derived by the Energetic Particle Detector (EPD) data on the Galileo mission, 1.1 * 10 26 atoms/s 

556 (Ipetal., 1998). 

557 Estimated energy spectra for IS, PSD, and ion backscattering (IBS) processes on Europa are shown 

558 in Figure A4 (Plainaki et al., 2010). Clearly, IS is far more productive, hence, SHEA dominates 

559 over releases from IBS (mainly H) and PSD at energies < 1 keV. 

560 The slightly lower incident ion fluxes and the similarity between Ganymede and Europa in surface 

561 composition, drives the conclusion that slightly less SHEA fluxes are expected at Ganymede, where 

562 the internal magnetic field is not able to shield the plasma (Kivelson et al. 2002). Callisto, on the 

563 other hand, is considerably out of Jupiter’s radiation belt; hence, the expected SHEA flux should be 



564 considerably lower in nominal conditions. A comparative detection of SHEA from these three 

565 Jovian moons would be of particular interest in the study of Jupiter system’s evolution. 

566 At Mercury, 1-keV solar wind H-ions release various types of atoms (like Na, Ca, K, Mg), and 

567 probably some molecules, too. Binding energies of these species with the surface of the planet are 

568 between 1 and few eVs. At the Galilean moons, 100-keV H, O and S ions of Jupiter's 

569 magnetospheric plasma, release mainly H 2 O. The sputtered particle energy distributions for 

570 molecular ices tend to have maxima at lower energies of about 0.05 eV (Boring et al. 1984, Haring 

571 et al., 1984). Simulations of ion sputtering show that at Europa and Mercury in the precipitation 

572 areas, die fluxes of the released particles differ at about one order of magnitude (10 9 particles/cm 2 /s 

573 at Europa (Plainaki et al., 2010) and 10 8 particles/cm 2 /s at Mercury (Mura et al., 2005)). However, 

574 according to Cassidy and Johnson (2005), in the non-ice regions of Europa, the regolith can 

575 significantly modify the relative populations of atmospheric species and their spatial distributions 

576 across the surface. Consequently, the sputtering yields should be reduced due to sticking of 

577 sputtered species to neighboring grains (Hapke, 1986; Johnson, 1989) and therefore lower fluxes of 

578 sputtered particles would be expected. 

579 Estimated escape fractions of sputtered particles from the different environments and the rough 

580 fraction of exposed bodies surfaces considered in this review are summarized in Table 1. 
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825 Table 1. Sputtered escape fraction from reference EB 


EB 

Exposed surface 
% 

Escaping sputtered fraction 
% 

Mercury 

25 

70 

Moon (out of Earth’s magnetosphere) 

50 

90 

NEO 

50 

100 

Europa 

100 

22 


826 

827 Figure 1 . Examples of ion-induced SHEA ejection. Incident energetic ions are in blue and atoms of 

828 the exposed surface are in red. Ejections shown are by (a) primary or first recoil, (b) knock-on by a 

829 backscattered ion without a cascade, (c) secondary recoil, (d) higher order recoil or cascade, (e) 

830 backscattered incident ion, and (f) backscattered recoil atom. Panels (e) and (f) are for surfaces 

831 having more than one element, e.g. 9 more massive atoms (in black). The angles shown in (c) define 

832 the directions of the incident ion, the recoil atom and the ejected SHEA. Panels (b) and (e) are 

833 indistinguishable externally, except in the energy of the SHEA. This figure is based on Figure 2.1 

834 of Sigmund (1981) and Figure 2.6 of Hofer (1991). 

835 Figure 2. Some measured and computed SHEA spectra, a) Computed energy distribution function 

836 fs t „ (Eq. (2)) of Na sputtered from 1-keV protons impacting on a simulated planetary- like mineral 

837 surface for different assumed Na binding energies, b) Computed /$„ of Na ejecta, assuming a 2 eV 

838 binding energy, for protons of different energies impacting on regolith-like simulated (again 

839 regolith means the porosity was account for), c) Comparison of high-energy part of Na ejecta 

840 distribution (dashed line) and SRIM simulations (solid line, for the assumed surface composition 

841 see Goettel, 1988) for Ej = 1 keV. d) SHEA energy spectra for Ar + of different E\ on W at zero 

842 incident and ejection angles (Goelhich et al. 2000). e) Ejection of Na from Na2S0 4 for impacting 

843 Ar + of 3.5 keV and with Eb - 0.27eV (Wiens et al. 1997). f) Sputtering of D 2 0 ice with - 30% S0 2 

844 bombarded by 5-keV Ar + . Energy profile of sputtered S0 2 (red dots) and D 2 0 (blue dots) 

845 molecules, normalized at 1 meV flux. The energy profiles are fit to two distributions of the form 

846 Eb/(E e +Eb) 2 . The fit shows for D0 2 (blue, lower curve), U ~ 0.048 eV for a fraction 0.32 of the 

847 molecules ejected with U ~ 0.003 3 eV for the remainder; for the S0 2 component (red, curve over the 




848 


dots), U ~ 0.043eV for a fraction 0.36 of the ejected molecules, with U ~ 0.0053 for the remainder 

849 (Johnson etal. 2010). 

850 Figure 3. Model of velocity distribution functions for TD (a), PSD (b) and IS (c) (adapted frn.: 

851 Killen et al., 2007). See also Fig.2 for SHEA spectra. Dashed lines correspond to the Na escape 

852 energy at Mercury equal to 2 eV, for reference. 

853 Figure Al. Top-left panel: Pseudo-color map of SFIEA flux from the Northern surface of Mercury 

854 due to proton sputtering (from Mura et al., 2005). The oval-shaped dotted line is the horizon as seen 

855 from the s/c (400 km above surface level); bottom-left panel: SHEA signal as seen from the s/c; 

856 bottom-right panel: ELENA FoV and count-rate (color) superimposed to the SHEA signal (gray); 

857 top-right panel: thanks to the s/c motion, the ground-track of ELENA data allows global imaging of 

858 the surface SHEA emission. 

859 Figure A2. Model of the equatorial lunar exosphere: (a) sodium density and its variation with solar 

860 zenith angle, %, and altitude for PSD, MIV, and IS; (b,c) sub-solar point profiles attributed to IS and 

861 MIV for a number of other abundant lunar constituents [From Sarontos et al., 2009], 

862 Figure A3. Simulated integral flux (log(particles m' 2 s' 1 )) of total sputtered particles from Cl 

863 chondrites NEO for impinging protons of energy ~1 keV. Axial symmetry is assumed; positive Y 

864 points to the Sun. (Plainaki et al., 2009). 

865 Figure A4. Intensity versus energy spectrum of the sputtered, back-scattered and PSD-ed neutrals at 

866 Europa (Plainaki et al., 2010) 
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